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Abstract. In this paper we demonstrate that radiation patterns could cause flow-like behaviour without 
any reference to hydrodynamic description. For that purpose we use a statistical ensemble of radiating 
dipoles, motivated by the investigation of the equivalent photon yield produced by decelerating charges. 

For the elliptic asymmetry factor, i> 2 , we find a reasonable agreement with experimental data. 
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Hydrodynamical simulations are widely used to de¬ 
scribe the early-time evolution of proton-nucleus and heavy- 
ion collisions, see Refs. mm and the references therein. 
However, the fact that hydrodynamics has a strong pre¬ 
dictive power does not imply that it is the only option to 
explain collective phenomena in such systems. There have 
been recent efforts to reproduce the flow patterns observed 
in RHIC and LHC using color scintillating antennas radi¬ 
ating gluons mm- Other authors utilized phenomeno¬ 
logical models of color-electric dipoles in order to account 
for angular correlations in high-energy processes mm. 
It is an ongoing debate though, whether a simple effec¬ 
tive model, lacking hydrodynamics, could catch the flow- 
like behaviour or not. Unfortunately, it is rather compli¬ 
cated to explain the collective properties using microscop¬ 
ical models as a starting point. 

In this letter, our goal is to demonstrate that the radia¬ 
tion originating from a dipole set-up is, in principle, able 
to match quantitatively the elliptic asymmetry factor V 2 , 
measured in heavy-ion experiments. To do so, we discuss 
the yield of massless particles produced by a decelerating 
point-like charge. Then we compute the flow coefficient i >2 
of a dipole composed of two, parallel displaced counter- 
decelerating charges. Motivated by the formula in Eq. d 6 j), 
we fit various experimental data. Finally, we conclude our 
analysis discussing several open issues and their relevance 
for further, more realistic description. 


Radiation produced by decelerating sources. According 
to electrodynamics, an accelerating point charge radiates. 
One can reinterpret this phenomenon quasi-classically as 
the emission of photons. It is straightforward to calculate 
the differential yield of emitted photons when the charge 


accelerates uniformly on a straight line m- 

W 2 k , 2 
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Here the azimuth angle if) is measured in the perpendic¬ 
ular plane, rj is the rapidity of the emitted photon in 
the laboratory frame, related to the detecting angle as 
tanh ?7 = cos 0, k± is the magnitude of the transverse mo¬ 
mentum of the photon. Parameters w\ and w 2 are the 
boosted initial and final velocities of the charge in the 
frame of the observer (Wi = 7 ivf), g is the magnitude of 
the co-moving acceleration (c = 1). It is noteworthy that 
Eq. JTJ in the k± —> 0 limit reproduces a bell-shaped ra¬ 
pidity distribution, similar to Landau’s hydrodynamical 
model,and also the plateau known from the Hwa-Bjorken 
scenario, depending whether the accelerating motion of 
the charge covers a short or large range in rapidity m- 
Speculating further, we assume that in the case of light 
particles produced in a heavy-ion collision a significant 
part of the yield comes from similar, deceleration induced 
radiation processes. 

It is also worthwhile to mention that a gauge field the¬ 
ory which describes the radiation phenomena on the mi¬ 
croscopic level, can be reformulated in the framework of 
hydrodynamics, as it was endeavoured in Refs. mm- 
Elliptic flow coefficient. We now turn to the analysis of 
the simplest structure. An elliptic asymmetry could stem 
from two decelerating point-like sources going into oppo¬ 
site directions on parallel paths. We calculate the emitted 
photon-equivalent radiation for the distance d between the 
two sources. The yield is calculated coherently summing 
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the contributing amplitudes: 

Y oc \Ai + A 2 1 2 = 

_ l^g^k-L cos(tp—il>) ^g-iffcx COs(tp—ll>) |2 

= I All 2 + \A 2 \ 2 + 2Re(A 1 A*e ik±dcosi *’-^). (2) 

Here p is the angle between the dipole axis and the detec¬ 
tor, and ip characterizes the event-plane. The elliptic flow 
coefficient we deal with here is defined by the azimuthal 
averaging of cos( 2 (y> — ip)) with respect to the normalized 
yield: 
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In the above formula J n is the Bessel function of the first 
kind. It is part of the Jacobi-Anger expansion, e l2:cos W = 
oo in Mz)e int , applied to e^ k± Integrat¬ 

ing with the factor cos( 2 (</j — ip)) and dividing by the ze¬ 
roth order coefficient delivers Eq. ©• Parametrizing the 
complex amplitudes as A\ = Ae l5 ° and A 2 = 'yAe l ^ 0+ ^ 
with real A , 7 and Sq, where <5 denotes the phase-shift 
and 7 is the ratio of |A 2 |/|Ai|, we obtain the simplified 
expression 

-J 2 (k±d)cosS 

V2 = TPT- T ,, n -7- ( 4 ) 

+ Jo(k±d) cos5 

Further we assume that A\ and A 2 are azimuthally sym¬ 
metric, i.e. independent of the difference ip — ip. 

In our picture v 2 depends on the phase-shift <5, the dipole 
size d and the strength asymmetry parameter 7 . We as¬ 
sume that event-by-event d and 7 might be well-determined, 
while d fluctuates. The decelerating sources are strongly 
affected by the medium surrounding them, therefore, they 
radiate differently, depending on how long the interaction 
holds up. We consider one source decelerating from a ve¬ 
locity near c to 0 and the other one from c to slightly above 
0, in the opposite direction. The relevant amplitudes then 
read: 
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as it follows from Eq. m with A = K n and L v being 
the modified Bessel functions of second kind and the mod¬ 
ified Struve function, respectively. The phase-shift factor, 
cos S = 1 can be evaluated numerically, the result 

is plotted in Fig. [ 1 ] It appears as a natural idea to average 
with respect to the phase difference variable, <5, whenever 
cos S oscillates fast as a function of w 2 (cf. Fig.©. 


Fig. 1: Phase-shift factor cosd as a function of the final 
velocity w 2 in cases A = 1.0 (dashed line) and A = 10.0 
(solid line). 


Fits to experimental data. The uniform averaging re¬ 
spect to the phase-shift angle can be carried out analyti¬ 
cally, resulting 
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Hereinafter we assume that the leading order contribution 
to the elliptic asymmetry comes from a dipole-like struc¬ 
ture discussed previously. We shall test our hypothesis on 
elliptic flow measurements in heavy-ion collisions of RHIC 
and LHC, where the fairly large number of dipoles ensures 
the validity of our working hypothesis, namely the uniform 
phase-shift of the sources. We introduce an additional fit 
parameter, called F. This geometrical form factor is in¬ 
dependent of the transverse momentum k±. Finally the 
formula we use to fit the experimental data scales the ex¬ 
pression © as: 


Mfit = F- {v 2 )(k±d,j). (7) 

Comparison of measured and calculated v 2 as a function 
of k± are depicted on Figs. © [J] and © For data resolved 
by centrality the fitted parameter, d , remains the same 
within 11% for photons and 19% for charged hadrons, cf. 
Table© Including more peripheral collisions, F saturates 
somewhat below 3 (see Fig©. Since F = 1 would mean 
that only the dipole term contributes to v 2 , cf. Eq. ©, 
that suggests the need for some further sources of elliptic 
flow, for example multipole contributions. In all cases 7 
turns out to be very close to one, showing that symmetric 
dipole sources may dominate the experimental findings. 
The v 2 values for charged pions on Fig. ©and hadrons on 
Fig. ©fit as well as for the emitted photons on Fig. © 
The typical value of the effective dipole size d is about 
0.06 fm from the hadronic fits and about 0.1 fm from 
photon data. This is rather small compared to the size 
of heavy-ion fireballs, it may hint to subhadronic sources 
of this part of radiation. We mention here that other au¬ 
thors pointed out the quark-level origin of the flow inde¬ 
pendently of our present analysis m- 
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Fig. 2: Inclusive photon elliptic flow measured by ALICE 
group of LHC in Pb-Pb collisions at ^/snn = 2.76 TeV 
for several centrality classes 1251 

At this point, the physical interpretation of the form fac¬ 
tor F is due. Being momentum independent, a simple ge¬ 
ometric cartoon of a heavy-ion collision can be suggested. 
Averaging the yield after Fourier-expansion of Eq. © with 
respect to a profile parametrized by polar coordinates we 
obtain 

Y = Y 0 (1) + Y 2 ( —\ cos 2ip = Y 0 + FY 2 cos 2 ip, (8) 

\ COS Zip / 

with ip being the polar angle measured from the reaction 
plane. In case of an ellipse, the factor in front of Y 2 is the 
anisotropy F = , with the half-axes A (in the reac¬ 

tion plane) and B (perpendicular to the reaction plane). 
There are several ways to attach an ellipse to the geometry 
of the collision. Let us consider two nuclei as circular disks 
(squeezed due to Lorentz contraction) with radius J?, dis¬ 
placed by impact parameter b between the centres. In case 
of approximating the intersection by an ellipse (A = R— |, 

B 2 = R 2 — - j-) and assuming subhadronic dipoles ordered 
perpendicular to the reaction plane we obtain F = 

A much larger ellipse including also the spectator area 
{A = R + |, B 2 = R 2 — ¥-) with dipoles ordered par¬ 
allel to the reaction plane delivers the same anisotropy 
F = 27 j. In fact experimentally F(b) turns out to be lin¬ 
ear in a wide range of impact parameter values, see Fig. 
[5l For fitting the anisotropy formula mentioned above we 
get R ~ 1.67fm. This is an effective size of the source of 
the dipole-like radiation. It is smaller than the typical size 
of a Pb-nucleus by a factor of four. This finding warns 
against a collective source extending in the whole media, 
but does not exclude hydrodynamic evolution. 

Discussion and conclusions. We briefly list recent lit¬ 
erature studies about how various stages of the heavy-ion 
collision could contribute to the azimuthal asymmetry of 
the flow in order to support the phenomenological picture 
we sketched above. We focus the role of dipole-like struc¬ 
tures revealing at the early-time stage of the heavy-ion 



k ± [GeV/c] 


Fig. 3: Elliptic flow of charged pions measured by STAR 
of RHIC in Au-Au collisions at yj snn = 62.4 GeV (solid 
line) [231 , an d by PHENIX of RHIC in Au-Au collisions 
at y/sNN = 200 GeV (dashed line) |2o] 



Fig. 4: Elliptic flow of charged hadrons measured by 
PHENIX in Au-Au collisions at y/s~NN = 200 GeV for 
several centrality classes m 

collisions. 

i.) Strong electromagnetic (EM) fields. In non-central col¬ 
lisions, the magnitude of the magnetic field due to the 
geometrical asymmetry of the system could reach ~ 5 m 2 
for a short time of 0.2fm/c [2HH2Sj • The pure EM-effect 
(caused by the coupling of charged quasi-particles and the 
EM-field) is, however, not significant at the level of global 
observables, as it is suggested by hadron string dynamics 
simulations mi or contributes to higher order asymme¬ 
tries only (quadrupole electric moment) [25]. Note that 
these simulations are based on transport models using 
quasi-particles and improved, but essentially perturbative 
cross sections. 

In Ref. m the authors use an order-of-magnitude estima¬ 
tion, leading order in perturbation theory, for the gluon- 
photon coupling in order to argue that the direct photon 
flow maybe affected at RHIC, but unlikely at LHC. 






































4 


T. S. Biro et al.: Elliptic flow due to radiation in heavy-ion collisions 


- 

•v''v V V V 

.-■Fm 




¥ ;y 


-M" 


photon form factor ■ 

M 

charged hadron form factor v 


linear fit to photons: 0.2938*b - 

\ 

\ 

linear fit to hadrons: 0.3026*b 


0 2 4 6 8 10 12 

Impact parameter, b [fm] 

Fig. 5: The form factor F - see Eq. 0 - versus the im¬ 
pact parameter b of the collision for inclusive photons and 
charged hadrons. The centrality-impact parameter rela¬ 
tionship is taken from m- 


ii.) QCD in magnetic field. Lattice Monte-Carlo simula¬ 
tions suggest that QCD at high temperatures is param¬ 
agnetic Ref. m- Therefore a ’squeezing’ of the plasma 
could occur, elongating it in the direction of external mag¬ 
netic field, which, in case of non-central collisions, points 
perpendicular to the reaction plane. Charge separation of 
quarks in the direction of the external magnetic field due 
to the fluctuation of the topological charge (known as the 
chiral magnetic effect, CME) can also contribute to the 
asymmetry of the plasma, as it is indicated by lattice re¬ 
sults |31] ■ These effects are not yet incorporated in quasi¬ 
particle simulations like in Refs. EEH22- 
in.) Radiation of non-Abelian plasma. The classical limit 
of non-Abelian fields generated by ultra-relativistic sources 
is analysed in Ref. ESI ■ R is shown that dipole-like struc¬ 
tures will emerge with the same geometric properties like 
their EM-versions. These could be important when initial 
state of the matter like the color glass condensate (CGC) 
is melting down and converting to QGP, while a consider¬ 
able amount of quark-antiquark pairs are produced. This 
happens probably when dipoles are smaller than 1 fm, ac¬ 
companied by fast oscillation of the sources. 

R is pointed out in Ref. [331 that the bremsstrahlung of 
quarks on the surface of QGP, pulled back by the confin¬ 
ing force could produce photon radiation in comparable 
amount to those produced in the plasma phase. An im¬ 
portant aspect of the issue of the azimuthal asymmetries 
is in what extent are those evolved on the microscopic 
level of the dynamics or caused by collective behaviour. 
There is an ongoing debate in the literature about the 
contributions of initial and final state asymmetries to the 
elliptic flow E3H22- Emphasizing only a few examples, 
it is observed that in proton-nucleon collisions the CGC- 
correlations could be directly visible in the measurable 
particle spectrum. Using classical Yang-Mills simulations 
for p-Pb collisions, in the first half fm/c after CGC was 
initiated, significant build-up of contributions to V 2 and 


V 3 was observed E3- These momentum space anisotropies 
are not correlated with the final state global asymmetries 
described as collective flow behaviour. Ri Pb-Pb collisions, 
the early-time contributions are relatively small, support¬ 
ing the role of collective effects. In this case the sources 
are uncorrelated, localized color field domains, resulting 
the gluon spectrum to be isotropic EH! • EM effects also 
could play a role in the final state. The directed flow of 
charged pions could be a result of a spectator induced 
splitting, as it is demonstrated in Refs. E3E2- 
It seems that dipole-like structures coupled to the initial 
geometric asymmetry of heavy-ion collisions are quite nat¬ 
ural in a wide scale of models concerning the early time 
evolution of the hot nuclear matter. We suggest that these 
domains could be the sources of intense gluon and/or pho¬ 
ton radiation having similar geometric properties to its 
EM counterparts. The orientation of these dipoles may 
be ordered by EM effects like the mentioned squeezing of 
the QCD plasma and CME, triggered by the early-time 
intense fields present in non-central events. Therefore the 
cumulative effect of small but not necessarily coherent ra¬ 
diators may affect the macroscopic observables, contribut¬ 
ing significantly to the azimuthally asymmetric compo¬ 
nent of the flow. 

In conclusion, we emphasize the necessity of exploring how 
microscopic causes can lead to macroscopic anisotropies. 
Especially in large systems, where collective (flow-like) ef¬ 
fects may take place, it is important to distinguish those 
from the individual subhadronic causes. 
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